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Liquid helium 4 had been the only bosonic superfluid available in experiments for a long time. 

This situation was changed in 1995, when a new superfluid was born with the realization of the 
Bose-Einstein condensation in ultracold atomic gases. The liquid helium 4 is strongly interacting 
and has no spin; there is almost no way to change its parameters, such as interaction strength and 
density. The new superfluid, Bose-Einstein condensate (BEC), offers various aspects of advantages 
over liquid helium. On the one hand, BEC is weakly interacting and has spin degrees of freedom. On 
the other hand, it is convenient to tune almost all the parameters of a BEC, for example, the kinetic 
energy by spin-orbit coupling, the density by the external potential, and the interaction by Feshbacli 
resonance. Great efforts have been devoted to studying these new aspects of superfluidity, and the 
results have greatly enriched our understanding of superfluidity. Here we review these developments 
by focusing on the stability and critical velocity of various superfluids. The BEC systems considered 
include a uniform superfluid in free space, a superfluid with its density periodically modulated, a 
superfluid with artificially engineered spin-orbit coupling, and a superfluid of pure spin current. Due 
to the weak interaction, these BEC systems can be well described by the mean field Gross-Pitaevskii 
theory and their superfluidity, in particular critical velocities, can be examined with Landau’s theory 
of superfluid. Experimental proposals to observe these new aspects of superfluidity are discussed. 

PACS numbers: 05.30.Jp, 03.75.Mn, 03.75.Kk, 71.70.Ej 


I. INTRODUCTION 

Superfluidity, as a remarkable macroscopic quantum 
phenomenon, was first discovered in the study of liq¬ 
uid helium 4 in the 1938 0,i- Although it is found 
theoretically that superfluidity is a general phenomenon 
for interacting boson systems Q , liquid helium had been 
the only bosonic superfluid available in experiments until 
1995. In this year, thanks to the advance of laser cooling 
of atoms, the Bose-Einstein condensation of dilute al¬ 
kali atomic gases was realized experimentally Q; a new 
superfluid, Bose-Einstein condensate (BEC), was born. 
This addition to the family of superfluids is highly non¬ 
trivial as BECs offer various aspects of advantages over 
liquid helium that can greatly enrich our understanding 
of superfluidity. 

Great deal of work has been done to explore the prop¬ 
erties of liquid helium as superfluid [1]. However, these 
efforts have been hindered by limitations of liquid helium. 
As a liquid, the superfluid helium is strongly-interacting 
system, which makes the theoretical description difficult. 
At the same time, no system parameters, such as density 
and interaction strength, can be tuned experimentally. 
And helium 4 has no spin degrees of freedom. 

BECs are strikingly different. Almost all the param¬ 
eters of a BEC can be controlled easily in experiments: 
its kinetic energy, density, and the interaction between 
atoms can all be tuned easily by engineering the atom- 
laser interaction, magnetic or optical traps, and the Fes- 
lrbach resonance @]. In addition, by choosing the atomic 
species and using optical trap to release the spin degrees 
of freedom, one can also realize various types of new 
superfluids, including the spinor superfluid [7| and the 


dipolar superfluid [§]. All these are impossible with liq¬ 
uid helium. Moreover, as most ultracold gases are dilute 
and weakly-interacting, controllable theoretical methods 
are available to study these superfluids in detail. 

In this review we mainly discuss three types of bosonic 
superfluids: superfluid with periodic density, superfluid 
with spin-orbit coupling, and superfluid of pure spin cur¬ 
rent. The focus is on the stability and critical velocities 
of various superfluids. For a uniform superfluid, when its 
speed exceeds a critical value, the system suffers Landau 
instability and superfluidity is lost. When the superfluid 
moves in a periodic potential, with large enough quasi¬ 
momentum, new mechanism of instability, i.e., dynamical 
instability comes in. This instability usually dominates 
the Landau instability as it occurs on a much faster time 
scale. The periodic density also brings another twist: the 
superfluidity can be tested in two different ways, which 
yield two different critical velocities. For a superfluid 
with spin-orbit coupling, a dramatic change is brought 
in, namely the breakdown of Galilean invariance. As a 
result, its critical velocity will depend on the reference 
frame. The stability of a pure spin current is also quite 
striking. We find that the pure spin current in general is 
not a super-flow. However, it can be stabilized to become 
a super-flow with quadratic Zeeman effect or spin-orbit 
coupling. Related experimental proposals are discussed. 

The paper is organized as follows. In Sec. UD we discuss 
briefly the basic concepts related to the understanding of 
superfluidity, including Landau’s theory of superfluidity, 
mean field Gross-Pitaevskii equation, and Bogoliubov ex¬ 
citations. These concepts are illustrated with the special 
case of a uniform superfluid. We then apply these gen¬ 
eral methods to study in detail the superfluid in periodic 
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potentials in Sec. um the superfluid with artificially en¬ 
gineered spin-orbit coupling in Sec. II VI and finally the 
superfluid of pure spin current in Sec. El In these three 
superfluids, special attention is paid to their excitations, 
stabilities, and critical velocities. We finally summarize 
in Sec. IVII 


its speed was smaller than the sound speed. This is in 
fact nothing but the Cerenkov radiation [To|, HH , where 
a charged particle radiates only when its speed exceeds 
the speed of light in the medium. 


II. BASIC CONCEPTS OF SUPERFLUIDITY 


A. Landau’s theory of superfluidity 


Superfluid is a special kind of fluid which does not 
suffer dissipation when flowing through a tube. It loses 
its superfluidity only when its speed exceeds a certain 
critical value. The superfluidity of liquid helium 4 was 
first explained by L. D. Landau [§]. He considered a 
superfluid moving inside a stationary tube with veloc¬ 
ity v. Since the system is invariant under the Galilean 
transformation, this scenario is equivalent to a stationary 
fluid inside a moving tube. If the elementary excitation 
in a stationary superfluid with momentum q has energy 
eo (q), then the energy of the same excitation in the back¬ 
ground of a moving fluid with v is e v (q) = eo(q)+ v 'R A 
fluid experiences friction only through emitting elemen¬ 
tary excitations, and it is a superfluid if these elementary 
excitations are energetically unfavorable. In other words, 
a superfluid satisfies the constraint e v (q) > 0. It readily 
leads to the well-known Landau’s criterion for superfluid, 


v <v c 



(1) 


Here v c is the critical velocity of the superfluid, which is 
determined by the smallest slope of the excitation spec¬ 
trum of a stationary superfluid. 

Another way of deriving the formula of critical velocity 
is from the point view of Cerenkov radiation. Consider a 
macroscopic impurity moving in the superfluid generates 
an excitation. According to the conservations of both 
momentum and energy, we should have 


m 0 Vi 

2 


m 0 v f + q, 


m 0 v 2 f 


+ eo(q), 


( 2 ) 

( 3 ) 


where too is the mass of the impurity, v, and v f are the 
initial and final velocities of the impurity, respectively. 
The above two conservations © and © can not be sat¬ 
isfied simultaneously when 


V « |Vj| « |v/| < v c 



( 4 ) 


The critical velocity v c here has the same expression as 
that obtained from Landau’s criterion. If the excitations 
are phonons, i.e., eo(q) = c|q|, then v c < c. This means 
that the impurity could not generate phonons in the su¬ 
perfluid and would not experience any viscosity when 


(a) (b) 

Superfluid V / \ Superfluid 



Superfluid 

V 


(c) 


FIG. 1: (a) A superfluid moves inside a stationary tube, (b) 
The superfluid is dragged by a tube moving at the speed of v. 
(c) An impurity moves at v in the superfluid. The two-way 
arrow indicates the equivalence between different scenarios. 

These two different ways of derivation are equivalent 
when the system has the Galilean invariance. By trans¬ 
forming to another reference frame illustrated in Fig. 
□2 b), the superfluid can be viewed as being dragged by 
a moving tube. We then replace the moving tube with 
a macroscopic impurity moving inside the superfluid as 
shown in Fig. ©c). Consequently, for systems with 
Galilean invariance, the critical velocity of a superfluid 
moving inside a tube without experiencing friction is just 
the critical velocity of an impurity moving in a superfluid 
without generating excitations. For systems without the 
Galilean invariance, such an equivalence is lost as we will 
see later with the spin-orbit coupled BEC. 

For ultracold bosonic gases, the low energy excitation 
is phonon, linear with respect to momentum q, thus the 
critical speed of this superfluid is just the sound speed. 
Taking into account the non-uniformity of the trapped 
gas, the critical velocity measured in experiments agrees 
well with the value predicted by Landau’s theory [121] . 
For liquid helium 4, as there is another kind of elemen¬ 
tary excitations called rotons, the critical speed of super- 
fluid helium 4 is largely determined by the roton excita¬ 
tion, much smaller than its sound speed. Nevertheless, 
the critical velocity measured in experiments is still one 
order of magnitude smaller than the value predicted by 
the theory [131 ] . So it is remarkable that Landau’s pre¬ 
diction of critical velocity was experimentally confirmed 
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with BEC almost six decades after its invention. 

One remark is warranted on Landau’s theory of su¬ 
perfluidity. Landau’s criterion © of critical speed does 
not apply for many superfluids. However, Landau’s ener¬ 
getic argument for superfluidity is very general and can 
be applied to all the cases considered in this review. We 
shall use this argument to determine the critical speeds 
of various superfluids. 


B. Gross-Pitaevskii equation and Bogoliubov 
excitation 

As liquid helium 4 is a strongly interacting system, 
the calculation of its excitation spectrum theoretically 
is challenging. However, for a dilute ultracold bosonic 
gases, convenient yet precise approximations can be made 
to determine its excitation in theory. Due to the dilute 
nature and short range interaction, the interaction be¬ 
tween atoms can be approximated by a contact interac¬ 
tion. Furthermore, for low energy scattering of bosonic 
atoms, the s-wave scattering channel dominates. So, as a 
good approximation, the actual complex interaction be¬ 
tween atoms is replaced by an effective s-wave contact 
interaction, i.e., U{ ri,r 2 ) = ^■Kh 2 a8(r\ — r 2 )/M, where a 
is the s-wave scattering length. At zero temperature, as¬ 
sume all the particles condense into the same orbit i/j { r ), 
then its evolution is governed by the mean field Gross- 
Pitaevskii (GP) equation [H[ , 

. ^CM) = -^^i’ + V{r)iP + c\ip\ 2 ip, (5) 

where V (r) is the trapping potential or other external 
potential, and c = 4-7T h 2 a/M is the interaction parameter 
or coupling constant. The GP energy functional reads 


into account the interaction between the condensed par¬ 
ticles. The first order approximation is to take into ac¬ 
count the interaction between the condensed particles 
and un-condensed particles. This is settled by the Bo¬ 
goliubov theory of elementary excitations The stan¬ 

dard Bogoliubov theory is to diagonalize the mean field 
Hamiltonian by the Bogoliubov transformation. Here we 
adopt an equivalent yet more convenient way to deal with 
this problem in inhomogeneous systems. For a stationary 
state ip(r,t) that satisfies the GP equation, we consider 
some small time dependent perturbations Sip(r,t) added 
to this stationary state. The perturbed state ’L = ip + Sip 
also satisfies the GP equation, 


ih 


9T(r, t) 
dt 


ft 2 V 2 
2 M 


T + P(r)T + c|T| 2 4'. 


(7) 


Expanding this equation to first order in the perturba¬ 
tion 8ip, one arrives at the equation of motion of the 
perturbation Sip(r,t), 


ih ~w = - ^ v2< ^ +y ( r W 

+2c\i/j\ 2 5if> + cip 2 Sip*, (8) 


and its complex conjugate partner. Assume the state 
before perturbation is a stationary state with the wave 
function ip{ r, t) = ^n (and write the perturba¬ 
tions as Sip(r,t) = \Jn(r)e~ lllt ^ h [u(r)e~ lult — v*(r)e lult ], 
then we arrive at a pair of equations that it(r) and c(r) 
satisfy, 

"($)-*($)■ (9) 

where the matrix A4 is given by 


S[W] 



~^\^\ 2 + v ( rM 2 + ^\ 4 


( 6 ) 

Note that here the wave function ip is normalized to the 
total particle number N. The approximation we make 
here is the mean field approximation, and the interac¬ 
tion between atoms is replaced by an effective mean field 
potential. The validity of the approximation usually de¬ 
pends on the condensed fraction. For weakly interact¬ 
ing dilute bosonic gases near zero temperature, the con¬ 
densed fraction can be more than 90 percent. In these 
situations, the GP equation works well. It has been 
used to calculate the collective excitations of trapped 
BEC, as well as vortex dynamics. Pretty good agree¬ 
ment is achieved between theory and experiment. With 
a stochastic term describing the effect of thermal atoms, 
the modified stochastic GP equation can also simulate 
BEC systems at finite temperature [n|. In addition, for 
atoms with spin, one can also derive a multicomponent 
GP equation to describe a spinor superfluid [lj|. 

The stationary state the GP equation describes is a 
zeroth order approximation in some sense. It only takes 


M = 


H 0 —cn(r) \ 
cn( r) -H 0 ) ’ 


( 10 ) 


with Ho = — ^yV 2 + P(r)— ^ + 2cn(r). By diagonalizing 
M, one obtains two sets of solutions, but only the solu¬ 
tion satisfying the constraint / dr (Im(i-)I 2 - |v(r)| 2 ) > 0 
represents physical excitations. Note that since the 
characteristic matrix A4 is not hermitian, its eigenval¬ 
ues are not necessarily real. If its eigenvalue e = hui 
has nonzero imaginary part, the perturbation Sip(r,t) = 
^/n(r)e _iMt / ? ‘[u(r)e _ “ t — v*(r)e lut ] will grow exponen¬ 
tially with time, signaling that the state before perturba¬ 
tion suffers dynamical instability. If its eigenvalue is real 
but negative, elementary excitations associated with the 
perturbations will be energetically favorable and super¬ 
fluidity is lost, which is called Landau instability. For the 
excitations of repulsive Bose gases without external po¬ 
tentials, the low energy excitations must be non-negative 
and gapless from the Hugenholtz-Pines theorem [181 ] . 

As a specific example, we apply the above formalism 
to the uniform Bose gas, namely, putting V (r) = 0. The 
wave function before perturbation is just ip = -yne ~ lfJ,t ^ h , 
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and the perturbation has this form u( r) = w q e lqr /\/V : , 
v(r) = r> q e zqr / y/V . Plugging these wave functions into 
the above equations, one immediately obtains the exci¬ 
tation spectrum for a uniform Bose gas, 

e i = \/ e q ( e ° + 2cn), (11) 

where = h 2 q 2 /2M is the single particle spectrum. At 
long wavelength or small momentum q , the excitation has 
the asymptotic form e q ~ qy/Wnc/M. This is nothing 
but phonon excitation, and \fnc/M is just the speed 
of sound. From Landau’s theory of critical velocity, we 
conclude that the critical velocity of a uniform superfluid 
Bose gas is just the sound speed. 

The method of mean field approximation and Bogoli- 
ubov transformation is very general, and applies in other 
more complicated situations. In the following discussion, 
we use this method to study three types of superfluids, 
superfluid in a periodic potential, superfluid with spin- 
orbit coupling and superfluid of pure spin current. 


III. PERIODIC SUPERFLUID 

It is hard to change the density of helium 4 as it is 
a liquid. In contrast, we can easily modulate the den¬ 
sity of a BEC which is a gas. When we put a BEC in 
an optical lattice, we obtain a superfluid whose density 
is periodically modulated. One can even further period¬ 
ically modulate the interatomic interaction of the BEC 
with optical Feshbach resonance |19]. Supersolid helium 
4 may be also regarded as a periodic superfluid as it can 
be viewed as some superfluid defects (most likely vacan¬ 
cies) flowing in a helium solid lattice 20]. In this section, 
we use a BEC in an optical lattice as an example to ex¬ 
amine the properties of a periodic superfluid. Compared 
to the uniform superfluid in free space, a new type of in¬ 
stability, i.e., the dynamical instability is found when the 
quasi-momentum k of the superfluid is larger than a crit¬ 
ical value. Usually the dynamical instability dominates 
the accompanying Landau instability as it happens on a 
much faster time scale [2l| . The presence of the periodic 
potential also brings along another critical velocity. 


A. Stability phase diagram 


Now we study the superfluidity of a BEC in a periodic 
potential (24 [23[, which is provided by the optical lattice 
in cold atom experiments. For simplicity, we consider a 
quasi-one-dimensional BEC, confined in a cigar-shaped 
trap. We treat the system with the mean field theory 
and obtain the grand-canonical GP Hamiltonian 


H = 



r 


Id 2 \ c. l4 

-2d^ +VC0SX U+2 W 


m 2 

( 12 ) 


where all the variables are scaled to be dimensionless with 
respect to a set of characteristic parameters of the sys¬ 
tem, the atomic mass M. the wave number kL of the 
laser light generating the optical lattice, and the aver¬ 
age density ng of the BEC. The chemical potential /i and 
the strength of the periodic potential v are in units of 
Ah 2 k 2 L /M, the wave function ip is in units of y/ng, x is in 
units of fci/2, and t is in units of M/Ah 2 k 2 L . The inter¬ 
action constant is given by c = 7rnoa s //c|, where a s > 0 
is the s-wave scattering length. 

For non-interacting case (c = 0), diagonalizing the 
Hamiltonian will give the standard Bloch waves and en¬ 
ergy bands. When the mean field interaction is turned 
on (c 0), in principle the Hamiltonian allows for 
other types of solutions which have no counterpart in 
the non-interacting case [U, [25[ . Here we focus on the 
solutions which still have the form of Bloch waves, i.e., 
ipk(x) = e lkx (pk{x), where cpk(x) has the same period with 
the optical lattice. 4>k(x ) can be found by extremizing the 
Hamiltonian above [231 ]. The solution found in this way 
should satisfy the stationary GP equation with periodic 
potential, 


1 d 2 

2 dx 2 


ip + v cos xip + c\ip\ ip = flip. (13) 


To determine the superfluidity of these Bloch states, 
we must consider elementary excitations around these 
Bloch states, and check whether the excitation energy is 
always positive. Positive excitation energy indicates that 
the Bloch state is a local energy minimum, and it is stable 
against small perturbations. Due to the periodicity of the 
Bloch wave, the perturbations can be decomposed into 
different decoupled modes labeled by q , 


S(pk(x,q) = u k {x,q)e lqx +v* k {x,q)e lqx , (14) 


where q ranges between —1/2 and 1/2 and the perturba¬ 
tion functions Uk and Vk are of periodicity of 27 t. 

Following the similar method in Sec. Ill Bl we linearize 
the GP equation above to obtain the Bogoliubov equation 
that Uk and Vk satisfy 


MM ( “‘ ) = e,(,) ( . (15) 


where 


M k (q) 


C(k + q) -c<p\ \ 

c&t 2 -C(~k + q) ) ’ 


(16) 


with 

\ f d \ 2 

C(k) = --( — +ikj + vcosx - fj, + 2c\(j) k \ 2 - (17) 


This eigenvalue equation has two sets of solutions, one 
corresponds to physical excitations, which is mostly 
phonon excitation, and the other can be called anti¬ 
phonon that is not physical. If the physical excitation 









v=0.01 


v=0.05 


v=0.10 


£fe(<?) is positive, the Bloch wave ipk is a local minimum, 
and the system will have superfluidity. Otherwise, the 
system suffers Landau instability or dynamical instabil¬ 
ity, depending on whether efc(g) is real negative or com¬ 
plex. 

In the case of free space v = 0, the Bloch state xpk 
becomes a plane wave e lkx . Then the operator A4^(q) 
becomes 


M k {q) = 


q 2 /2 + kq + c — c 

c —q 2 /2 + kq~c 


(18) 


and we recover the excitations in the uniform case 


£±(q) =kq±^cq 2 + g 4 /4. (19) 


One immediately sees that the excitation energy is always 
real, which means that the BEC flows in free space are 
always dynamically stable. 

When there is periodic potential, the situation is dra¬ 
matically different, where the excitation energy can have 
imaginary part, signaling the dynamical instability of the 
system. By numerically solving the Bogoliubov equation 
above, we show the stability phase diagrams for BEC 
Bloch waves in the panels of Fig. [2j where different val¬ 
ues of v and c are considered. The results have reflection 
symmetry in k and q , so we only show the parameter 
region, 0 < k < 1/2 and 0 < q < 1/2. In the shaded 
area (light or dark) of each panel of Fig. [2J the exci¬ 
tation energy is negative, and the corresponding Bloch 
states ipk are saddle points. For those values of k out¬ 
side the shaded area, the Bloch states are local energy 
minima and represent superfluids. The superfluid region 
expands with increasing atomic interaction c, and occu¬ 
pies the entire Brillouin zone for sufficiently large c. On 
the other hand, the lattice potential strength v does not 
affect the superfluid region very much as we see in each 
row. The phase boundaries for v <C 1 are well reproduced 
from the analytical expression k = \Jq 2 1 4 + c for v = 0, 
which is plotted as triangles in the first column. 

If ek(q) is complex, the system suffers dynamical in¬ 
stability, which is shown by the dark-shaded areas in 
Fig. [2] The dynamical instability is the result of the 
resonance coupling between a phonon mode and an anti¬ 
phonon mode by first-order Bragg scattering. The matrix 
A4k{q) is real in the momentum representation, meaning 
that its complex eigenvalue can appear only in conjugate 
pairs and they must come from a pair of real eigenvalues 
that are degenerate prior to the coupling. Degeneracies 
or resonances within the phonon spectrum or within the 
anti-phonon spectrum do not give rise to dynamical in¬ 
stability; they only generate gaps in the spectra. Based 
on this general conclusion, we consider two special cases, 
which allow for simple explanations of the onset of dy¬ 
namical instability. 

One case is the weak periodic potential limit » « 1, 
where we can approximate the boundary with the free 
space case. This case corresponds to the first column of 
Fig. [2] In this limit, we can approximate the phonon 



FIG. 2: Stability phase diagrams of BEC Bloch states in op¬ 
tical lattices, k is the wave number of BEC Bloch waves; 
q denotes the wave number of perturbation modes. In the 
shaded (light or dark) area, the perturbation mode has neg¬ 
ative excitation energy; in the dark shaded area, the mode 
grows or decays exponentially in time. The triangles in (a.l- 
a.4) represent the boundary, q 2 /A + c = fc 2 , of saddle point 
regions at v = 0. The solid dots in the first column are from 
the analytical results of Eg. (1201) . The circles in (b.l) and (c.l) 
are based on the analytical expression m■ The dashed lines 
indicate the most unstable modes for each Bloch wave tpk- 


spectrum and the anti-phonon spectrum with the ones 
given in Eq. m- By equating them, e+(g — 1) = e_(g), 
for the degeneracy, we find that the dynamical instability 
should occur on the following curves 

k = \/q 2 c + g 4 /4 + i/(g - l) 2 c + (q - l) 4 /4 . (20) 

These curves are plotted as solid dots in Fig. [21 and they 
fall right in the middle of the thin dark-shaded areas. To 
some extent, one can regard these thin dark-shaded areas 
as broadening of the curves (l20l) . It is noted in Ref. IHj 
that the relation CHI) is also the result of e + (q — 1) + 
e_|_(— q) = 0, which involves only the physical phonons. 
Therefore, the physical meaning of Eq. (l20l) is that one 
can excite a pair of phonons with total energy zero and 
with total momentum equal to a reciprocal wave number 
of the lattice. 

The other case, c <C v, is shown in the first row of 
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Fig. [2] The open circles along the left edges of these 
dark-shaded areas are given by 

Ei(k + q)~ E\{k) = Ei(k) - E^k - q), (21) 

where E\(k) is the lowest Bloch band of 

1 <9 2 

H ° = + vcos ^- ( 22 ) 

In this linear periodic system, the excitation spectrum 
(phonon or anti-phonon) just corresponds to transitions 
from the Bloch states of energy E\{k) to other Bloch 
states of energy E n {k + q), or vise versa. The above 
equation is just the resonance condition between such 
excitations in the lowest band (n = 1). Alternatively, we 
can write the resonance condition as 

£i(fc) + E^k) = E 1 {k + q) + E\(k - q) . (23) 

So, this condition may be viewed as the energy and mo¬ 
mentum conservation for two particles interacting and 
decaying into two different Bloch states E\{k + q) and 
Ei {k — q). This is the same physical picture behind Eq. 

®. 

One common feature of all the diagrams in Fig. [2] is 
that there are two critical Bloch wave numbers, k t and 
kd. Beyond kt the Bloch waves il>k suffer the Landau in¬ 
stability; beyond kd the Bloch waves i pk are dynamically 
unstable. The onset of instability at kd always corre¬ 
sponds to q = 1/2. In other words, if we drive the Bloch 
state il>k from fc = 0tofc = l/2 the first unstable mode 
appearing is always q = ±1/2, which represents period 
doubling. Only for k > kd can longer wavelength instabil¬ 
ities occur. The growth of these unstable modes drives 
the system far away from the Bloch state and sponta¬ 
neously breaks the translational symmetry of the system. 


B. Two critical velocities 

Besides inducing the dynamical instability, the pres¬ 
ence of the optical lattice has also non-trivial conse¬ 
quences on the concept of critical velocity. In contrast to 
the homogeneous superfluid which has only one critical 
velocity, there are two distinct critical velocities for a pe¬ 
riodic superfluid |27[. The first one, which we call inside 
critical velocity, is for an impurity to move frictionlessly 
in the periodic superfluid system (Fig. [3])a)) ; the second, 
which is called trawler critical velocity, is the largest ve¬ 
locity of the lattice for the superfluidity to maintain (Fig. 
[3Kb)). We illustrate these two critical velocities with a 
BEC in a one-dimensional optical lattice. 

The presence of the optical lattice plays a decisive role 
in the appearance of the two critical velocities: two very 
different situations can arise. The first situation is de¬ 
scribed in Fig. m a )> where one macroscopic impurity 
moves inside the superfluid. The key feature in this situ¬ 
ation is that there is no relative motion between the su¬ 
perfluid and the lattice. The other situation is illustrated 



FIG. 3: (a) A macroscopic obstacle moves with a velocity of 
v inside a superfluid residing in a periodic potential. The 
superfluid and periodic potential are “locked” together and 
there is no relative motion between them, (b) The lattice 
where a superfluid resides is slowly accelerated to a velocity 
of v. 


in Fig. EKb), where the lattice is slowly accelerated to a 
given velocity and there is a relative motion between the 
superfluid and the lattice. For these two different situa¬ 
tions, naturally arise two critical velocities. 

In the first scenario, we consider a moving impurity 
that generates an excitation with momentum q and en¬ 
ergy e 0 (q) in the BEC. According to the conservations of 
both momentum and energy, we should have 


m 0 Vi 

w 0 vf 

2 


mov f + q + nh.G , 


TOoV 2 


eo(q) ■ 


(24) 

(25) 


where mg is the mass of the impurity, Vi and v f are the 
initial and final velocities of the impurity, respectively, 
G is the reciprocal vector, and eo(q) is the excitation of 
the BEC at the lowest Bloch state k = 0. Note that in 
contrast to the conservation of momentum in free space, 
here the total momentum of the impurity and the ex¬ 
citation is not exactly conserved due to the presence of 
an optical lattice, and the momenta differing by integer 
multiples of reciprocal lattice vector are equivalent. For 
phonon excitations, i.e., eo(q) = c|q|, these two conser¬ 
vations can always be satisfied simultaneously no matter 
how small the velocity of the impurity is. In other words, 
the critical velocity for this scenario is exactly zero. 

In the other scenario, there is relative motion between 
the superfluid and the optical lattice, and the superfluid 
no longer resides in the k = 0 point of the Brillouin 
zone. We should examine the stability of Bloch waves 
with nonzero fc, which is discussed in detail in the previ¬ 
ous subsection. The critical wave number k t mentioned 
above corresponds precisely to the trawler critical veloc¬ 
ity Vt here. As kt is not zero, i>t is not zero. 

Both critical velocities can be measured with BECs 
in optical lattices. The inside critical velocity Uj can 
be measured with the same experimental setting as in 
Ref. & where the superfluidity of a BEC was stud¬ 
ied by moving a blue-detuned laser inside the BEC. For 
the trawler critical velocity Vt , one can repeat the exper- 






7 


iment in Ref. (2lj . where a BEC is loaded in a moving 
optical lattice. One only needs to shift his attention from 
dynamical instability to superfluidity. The potential dif¬ 
ficulty lies in that the Landau instability occurs over a 
much larger time scale, which may be beyond the life 
time of a BEC 0- 

IV. SUPERFLUIDITY WITH SPIN ORBIT 
COUPLING 

The intrinsic spin-orbit coupling (SOC) of electrons 
plays a crucial role in many exotic materials, such as 
topological insulators [29|. In spintronics [30], its pres¬ 
ence enables us to manipulate the spin of electrons by 
means of exerting electric field instead of magnetic field, 
which is much easier to implement for industrial applica¬ 
tions. However, as a relativistic effect, the intrinsic SOC 
does not exist or is very weak for bosons in nature. With 
the method of engineering atom-laser interaction, an ar¬ 
tificial SOC has been realized for ultracold bosonic gases 
in PS. A great deal of effort has been devoted to 
study many interesting properties of spin-orbit coupled 
BECs PP]. 

A dramatic change that the SOC brings to the concept 
of superfluidity is the breakdown of Landau’s criterion of 
critical velocity © and the appearance of two different 
critical velocities. Laudau’s criterion of critical velocity 
© is based on the Galilean invariance. It is apparent 
to many that the scenario where a superfluid flows in¬ 
side a motionless tube is equivalent to the other scenario 
where a superfluid at rest is dragged by a moving tube. 
If the flowing superfluid loses its superfluidity when its 
speed exceeds a critical speed v c , then the superfluid in 
the other scenario will be dragged into motion by a tube 
moving faster than v c . However, this equivalence is based 
on that the superfluid is invariant under the Galilean 
transformation. As SOC breaks the Galilean invariance 
of the system [48|, we find that the two scenarios men¬ 
tioned above are no longer equivalent as shown in Fig. [I] 
the critical speed for scenario (a) is different from the one 
for scenario (b). For easy reference, the critical speed for 
(a) is hereafter called the critical flowing speed and the 
one for (b) the critical dragging speed. 

For ultra-cold atomic gases, the breakdown of the 
Galilean invariance at the presence of SOC can be un¬ 
derstood both theoretically and experimentally. 

From the theoretical point of view, we show in de¬ 
tail how a system with SOC changes under the Galilean 
transformation. We adopt the formalism in Ref. [48| . 
The operator for the Galilean transformation is 


G(v, t) = exp [zv • (mr — pt)/h] , (26) 

which satisfies the definition 

G^(v,i)rG(v,f) = r — vi, (27) 

G t (v,t)pG(v,t)= p — mv, (28) 

G*(v, t)<xG(v, i) = er. (29) 


(a) 


(b) 


SOC superfluid 


<7^ 


SOC superfluid 


V d 



11 s 

•JJ' 



(0 


FIG. 4: (a) A superfluid with SOC moves while the tube is at 
rest, (b) The superfluid is dragged by a tube moving at the 
speed of v. (c) An impurity moves at v in the SOC superfluid. 
The reference frame is the lab. The two-way arrow indicates 
the equivalence between different scenarios and the arrow with 
a bar indicates the non-equivalence. 


A system is invariant under the Galilean transformation 
if the following equation is satisfied (Ref. 0), 


G+(v,«) 


‘ h ^- H 


G(v,t) = 


d 

ih dt~ H 


(30) 


The above condition is clearly satisfied by a Hamiltonian 
without SOC, e.g., H = p 2 /2m. However, for a Hamil¬ 
tonian with SOC, e.g., H soc = p 2 /2m + 7 (a x p y - a v p x ), 
it is easy to check that 


G t (v, t) 


(h 777 H soc 

at 


G(v,t) = 


4 - K. 


. (31) 


where H' soc = H soc + my (a x v y — a y v x ). Clearly there 
is an additional term dependent on the velocity of the 
reference frame. This new term can be regarded as an 
effective Zeeman effect and can not be gauged away; the 
Galilean invariance of the system is thus lost. 

In the experiments of ultra-cold atomic gases, the SOC 
is created by two Raman beams that couple two hyperfine 
states of the atom. Since the Galilean transformation 
only boosts the BEC, not including the laser setup as a 
whole, the moving BEC will experience a different laser 
field due to the Doppler effect, resulting a loss of the 
Galilean invariance. 


A. Bogoliubov excitations and definition of critical 
velocities 

We use the method introduced in Sec. IIIBI to study 
the superfluidity of a BEC with SOC by computing its el¬ 
ementary excitations Rol . In experiments, only the equal 





















combination of Rashba and Dresselhaus coupling is real¬ 
ized. Here we use the Rashba coupling as an example as 
the main conclusion does not rely on the details of the 
SOC type. 

We calculate how the elementary excitations change 
with the flow speed and manage to derive from these 
excitations the critical speeds for the two different sce¬ 
narios shown in Fig. 0Ja,b). We find that there are two 
branches of elementary excitations for a BEC with SOC: 
the lower branch is phonon-like at long wavelengths and 
the upper branch is generally gapped. Careful analysis of 
these excitations indicates that the critical flowing veloc¬ 
ity for a BEC with SOC (Fig. |4|a)) is non-zero while the 
critical dragging speed is zero (Fig. QJb)). This shows 
that the critical velocity depends on the reference frame 
for a BEC with SOC and, probably, for any superfluid 
that has no Galilean invariance. 

Specifically, we consider a BEC with pseudo-spin 1/2 
and the Rashba SOC. The system can be described by 
the Hamiltonian [3(| 39;, 50, .51] 


H = f dr { E 2 ^(-w +y(r) )^ 

+7 bPl (ifix + Py)l !)2 + tp2 {-iPx + Py)^l] 

+Y (IV’il 4 + IV^I 4 ) + c 2 |V’il 2 IV’ 2 | 2 | , (32) 

where 7 is the SOC constant, Ci and C 2 are interaction 
strengths between the same and different pseudo-spin 
states, respectively. For simplicity and easy comparison 
with previous theory, we focus on the homogeneous case 
F(r) =0 despite that the BEC usually resides in a har¬ 
monic trap in experiments. Besides, we limit ourselves 
to the case Ci > C 2 , namely, in the plane wave phase. In 
the following discussion, for simplicity, we set H = M = 1 
and ignore the non-essential z direction, treating the sys¬ 
tem as two-dimensional. We also assume the BEC moves 
in the y direction, and the critical velocity is found to be 
not influenced by the excitation in the z direction. 

The GP equation obtained from the Hamiltonian (l32l) 
has plane wave solutions 


0k = 



1 ( e idk 
- 1 


ikr—i^(k)£ 


(33) 


where tan 8^ = k x /k y , k) = |k| 2 / 2 — 7 |k| + (C'i + C , 2 )/ 2 . 
The solution 0k is the ground state of the system when 
|k| = 7 . There are another set of plane wave solutions, 
which have higher energies and are not relevant to our 
discussion. 

We study first the scenario depicted in Fig. |U]a), where 
the BEC flows with a given velocity. Since the system 
is not invariant under the Galilean transformation, we 
cannot use Laudau’s argument to find the excitations for 
the flowing BEC from the excitation of a stationary BEC. 
We have to compute the excitations directly. This can 
be done by computing the elementary excitations of the 


state 0k with the Bogoliubov equation for different values 

of k. 

Without loss of generality, we choose k = ky with k > 
0. Following the standard procedure of linearizing the 
GP equation [22|, |23[, we have the following Bogoliubov 
equation 


[ U1 \ 

( u l\ 

U2 \=e 

U 2 

V! 

Vi 

\ t ’2 / 

\v 2 J 


(34) 


where 


Ht 

b \2 

-hCi 


&21 


-2C2 

-he, 

hci 

hc-2 

u k 

^34 


iCi 

^43 

u k 


(35) 


with = ± gg+(9v ± fc ) 2 ± a, A = % - 4 + 7*, bi 2 = 
- 7 (iQx + q v + k) + fy, 621 = 7 (*<?x - q y - k) + , 634 = 

7 (iqx ~ Q v + k) - and 643 = -7 + q y - k) - &. 
As usual, there are two groups of eigenvalues and only 
the ones whose corresponding eigenvectors satisfy |iq | 2 — 
|ui | 2 = 1 (i = 1 , 2 ) are physical. 

In general there are no simple analytical results. We 
have numerically diagonalized A4 to obtain the elemen¬ 
tary excitations. We find that part of the excitations 
are imaginary for BEC flows with |k| < 7 . This means 
that all the flows with |k| < 7 are dynamically unstable 
and therefore do not have superfluidity. For other flows 
with k > 7 , the excitations are always real and they are 
plotted in Fig. 0 One immediately notices that the ex¬ 
citations have two branches, which contact each other at 
a single point. Closer examination shows that the upper 
branch is gapped in most of the cases while the lower 
branch has phonon-like spectrum at large wavelength. 
These features are more apparent in Fig. [Gl where only 
the excitations along the x axis and y axis are plotted. 

In Fig. [Sjc) and Fig. [Htc2), we notice that part of the 
excitations in the upper branch are negative, indicating 
that the underlying BEC flow is thermodynamically un¬ 
stable and has no superfluidity. In fact, our numerical 
computation shows that there exists a critical value k c : 
when k > k c either part of the upper branch of excita¬ 
tions or part of the lower branch or both become nega¬ 
tive. This means that the flows described by the plane 
wave solution 0k,- with |k| > k c suffer Landau instabil¬ 
ity and have no superfluidity. Combined with the fact 
that the flows with |k| < 7 are dynamically unstable, 
we can conclude that only the flows with 7 < |k| < k c 
have superfluidity. The corresponding critical speed is 
v c = k c — 7 . We have plotted how the critical flowing 
velocity v c varies with the SOC parameter 7 in Fig. [7] 

We turn to another reference frame illustrated in Fig. 
|4jb), where the BEC can be viewed as being dragged by 
a moving tube. To simplify the discussion, we replace 
the moving tube with a macroscopic impurity moving in¬ 
side the BEC as shown in Fig. @|c). Correspondingly, 
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FIG. 6 : Excitations along the x axis (the first row) and y 
axis (the second row) at different values of k. (al,a 2 ) k — 1 ; 
(bl,b 2 ) k = 3; (cl,c 2 ) k = 4. Ci = 10 , C 2 = 4, 7 = 1 . 


(a) 


(b) 




FIG. 7: Critical flowing velocity v c of a BEC as a function of 
the SOC parameter 7 . (a) Ci = 11, C 2 = 3; (b) Ci = 14, 
C 2 = 6 . 


the question “whether the BEC will be dragged along by 
the moving tube?” is replaced by an equivalent question 
“whether the impurity will experience any viscosity?”. 
Suppose that the moving impurity generates an excita¬ 
tion in the BEC. According to the conservations of both 


momentum and energy, we should have 


moVi = mov/ + q, 


w 0 v.f 

2 


m 0 -v 2 f 


+ e 0 (q), 


(36) 

(37) 


as those for a BEC without SOC, where £o(q) is the ex¬ 
citation of the BEC at k = 7 . The critical dragging 
velocity derived from Eqs. (1361371) is given by 


Vc 


ep(q) 

|q| 


(38) 


If the excitations were purely phonons, i.e., eo(q) = c|q|, 
these two conservations would not be satisfied simulta¬ 
neously when v « | | ss | vy | < c. This means that the 

impurity could not generate phonons in the superfluid 
and would not experience any viscosity when its speed 
was smaller than the sound speed. Unfortunately, for 
our BEC system, the elementary excitations eo(q) are 
not purely phonons, as will be shown below. 

When 7 /O, the excitations eo(q) also share two 
branches. Along the x axis, these two branches are 



(39) 

where Si = 2y 4 + 7 2 (Ci - C2), S 2 = 2 y 2 + \C\, t\ = sf, 

t 2 = 2 siS 2 , and t .3 = 2si + (y 2 + C 2 / 2 ) 2 ■ Along the y 
axis, the excitations of the ground state are 

\ /dd+C 2 9 , q$ 

e o( 9 y) = y—2— q y + 4 ’ ( 4 °) 

4 (%) = 2 7% + 2s ! + (^ s 2 - y ) q y + ' J - ( 41 ) 

When 7 > 0, the upper branch (q x ) is always parabolic 
at small q x with a gap \/2si. When expanded to the 
second order of q x , the lower branch has the following 
form 

£o (<?*) » q 2 x \J Cl ^ 2 ■ ( 42 ) 
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This shows that e/"(9®) is parabolic at long wavelengths 
instead of linear as usually expected for a boson system. 
This agrees with the results in Ref. [4l|. This parabolic 
excitation has a far-reaching consequence: according to 
Eq. (l38l) . the critical dragging velocity v c is zero, very 
different from the critical flowing velocity for a BEC mov¬ 
ing in a tube. This shows that the critical velocity for a 
BEC with SOC is not independent of the reference frame, 
in stark contrast with a homogeneous superfluid without 
SOC. This result of course has the root in the fact that 
the BEC described by the SOC Hamiltonian (15^1) is not 
invariant under the Galilean transformation [481 1 - 

We have also investigated the superfluidity with the 
general form of SOC, which is a mixture of Rashba and 
Dresselhaus coupling. Mathematically, this SOC term 
has the form aa x p y — /3a y p x . The essential physics is 
the same: the critical flowing speed is different from the 
critical dragging speed, and therefore the critical velocity 
depends on the choice of the reference frame. However, 
the details do differ when a ^ /3. The critical dragging 
speed is no longer zero. Without loss of generality, we 
let a > p. The slope of the excitation spectrum for the 
ground state along the y axis is 


v y = ^\/2a 2 (Ci - C 2 + 2a 2 ) + 2a 2 + Cl g ° 2 
and the slope along the x axis, 


— 2 a , 
(43) 


v x 



C\ + C"2 
2 


(44) 


The critical dragging velocity is the smaller one of the 
above two slopes, both of which are nonzero. 


B. Experimental observation 

Spin-orbit couple_d BECs have been realized by many 
different groups [dll - lddl] through coupling ultracold 
8 'Rb atoms with laser fields. The strength of the SOC in 
the experiments can be tuned by changing the directions 
of the lasers [3ll - l33j or through the fast modulation of 
the laser intensities [52|. The interaction between atoms 
can be adjusted by varying the confinement potential, the 
atom number or through the Feshbach resonance Q. For 
the scenario in Fig. Q|b), one can use a blue-detuned laser 
to mimic the impurity for the measurement of the critical 
dragging speed similar to the experiment in Ref. fib? . For 
the scenario in Fig. Ufa), there are two possible exper¬ 
imental setups for measuring the critical flowing speed. 
In the first one, one generates a dipole oscillation similar 
to the experiment in Ref. |33] but with a blue detuned 
laser inserted in the middle of the trap. The second one 
is more complicated: At first, one can generate a moving 
BEC with a gravitomagnetic trap [54|. One then uses 
Bragg spectroscopy [53, |5g to measure the excitations 
of the moving BEC, from which the superfluidity can be 


inferred. For the typical atomic density of 10 14 ~ 10 15 
cnr 3 achievable in current experiments (53j, and the ex¬ 
perimental setup in Ref. JUJ, the critical flowing velocity 
is 0.2 ~ 0.6 mm/s, while the critical dragging velocity 
is still very small, about 10 -3 ~ 10 -2 mm/s. To further 
magnify the difference between the two critical velocities, 
one may use the Feshbach resonance to tune the s-wave 
scattering length. 


V. SUPERFLUIDITY OF SPIN CURRENT 

A neutral boson can carry both mass and spin; it thus 
can carry both mass current and spin current. However, 
when a boson system is said to be a superfluid, it tra¬ 
ditionally refers only to its mass current. The historical 
reason is that the first superfluid discovered in experi¬ 
ment is the spinless helium 4 which carries only mass 
current. For a boson with spin, say, a spin-1 boson, we 
can in fact have a pure spin current, a spin current with 
no mass current. This pure spin current can be gener¬ 
ated by putting an unpolarized spin-1 boson system in 
a magnetic field with a small gradient. Can such a pure 
spin current be a super-flow? In this section we try to 
address this issue for both planar and circular pure spin 
currents by focusing on an unpolarized spin-1 Bose gas. 

The stability of a pure spin current in an unpolarized 
spin-1 Bose gas was studied in [57} ■ It was found that 
such a current is generally unstable and is not a super- 
flow. We have recently found that the pure spin current 
can be stabilized to become a super-flow |58[: (*) for a 
planar flow, it can be stabilized by quadratic Zeeman ef¬ 
fect; (it) for a circular flow, it can be stabilized with SOC. 
We shall discuss these results and related experimental 
schemes in detail in the next subsections. 

There has been lots of study on the counterflow in a 
two-species BEC [H3, lHHiH , which appears very similar 
to the pure spin current discussed here. It was found 
that there is a critical relative speed beyond which the 
counterflow state loses its superfluidity and becomes un¬ 
stable HHszl. We emphasize that the counterflow in a 
two-species BEC is not a pure spin current for two rea¬ 
sons. Firstly, although the two species may be regarded 
as two components of a psudo-spin 1/2, they do not have 
SU(2) symmetry. Secondly, it is hard to prepare experi¬ 
mentally a BEC with exactly equal numbers of bosons in 
the two species to create a counterflow with no mass cur¬ 
rent. There is also interesting work addressing the issue 
of spin superfluidity in other situations j70l - 172| . 


A. Planar flow 


We consider a spin-1 BEC in free space. The mean- 
field wave function of such a spin-1 BEC satisfies the 
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following GP equation [73}, 

q S, 2 V 2 1 v 

ih—if) m = — ^-V , m + co7'V , m + C 2 ^ s-S mn i/j n , (45) 

~ n=— 1 


where (to = 1 , 0 ,— 1 ) are the components of the 
macroscopic wave function, p = Y^ m =-i |'*/’m | 2 is the 
total density, s» = J 2 mn' t l ) m(Si)mnipn is the spin density 
vector and S = (S x , S y , S z ) is the spin operator vector 
with Si ( i = x, y, z) being the three Pauli matrices in 
the spin-1 representation. The collisional interactions in¬ 
clude a spin-independent part Cq = 47 rfi. 2 (ao + 2 a 2 )/ 3 M 
and a spin-dependent part C 2 = 47 rfi, 2 (a 2 — ao)/3 M, with 
a/(/ = 0 , 2 ) being the s-wave scattering length for spin -1 
atoms in the symmetric channel of total spin /. 

We consider a spin current state of the above GP equa¬ 
tion with the form 


= 



(46) 


where n is the density of the uniform BEC. The require¬ 
ment of equal chemical potential leads to |ki| = | k 2 1 - In 
the case where ki = ~k 2 , this state carries a pure spin 
current: the total mass current is zero as it has equal 
mass counterflow while the spin current is nonzero. 

To determine whether the state (l46l) represents a su¬ 
perfluid, we need compute its Bogoliubov excitation spec¬ 
trum, also using the method introduced in Sec. Ill Bl 
It is instructive to first consider the special case when 
there is no counterflow, i.e., ki = k 2 = 0. The exci¬ 


tation spectra are found to be e° = \j^C 2 ne q + e 2 and 
ef 1 = yj 2 c 0 ne q + e 2 , e ^ 1 = ^ 2 c 2 ne q + e 2 , respectively, 


with e q = h 2 q 2 /2M. So for antiferromagnetic interac¬ 
tion (co > 0 , C 2 > 0 ), all branches of the spectra are real 
and there is a double degeneracy in one branch of the 
spectra. The phonon excitations give two sound veloci¬ 
ties, \JnCi/M (i = 0 , 2 ), corresponding to the speeds of 
density wave and spin wave, respectively. However, the 
existence of phonon excitation does not mean that the 
pure spin current (ki = k 2 7 ^ 0 ) is a super-flow as we can 
not obtain the current with ki = k 2 7 ^ 0 from the state 
with ki = k 2 = 0 by a Galilean transformation. 

For the counterflow state with ki = —k 2 7 ^ 0, the sta¬ 
bility has been studied in Ref. [57|. It is found that, for 
the antiferromagnetic interaction case (co > 0 ,C 2 > 0 ), 
the excitation spectrum of the m = 0 component always 
has nonzero imaginary part in the long wavelength limit 
as long as there is counterflow between the two compo¬ 
nents, and the imaginary excitations in the to = 1,-1 
components only appear for a large enough relative ve¬ 
locity v\ = 2 yJn,C 2 /M. For the ferromagnetic interac¬ 
tion case (co > 0 , C 2 < 0 ), both excitation spectra of the 
to = 0 and to = 1,-1 components have nonzero imagi¬ 
nary parts for any relative velocity. This means that the 
pure spin current cannot be stable in any cases. 


For the general non-collinear case (k = kl ^ k2 7 ^ 0) 
and antiferromagnetic interaction, the excitation spec¬ 
trum for the to = 0 component is found to be 


e° = \j(e q + ^ (|k| 2 - |ki| 2 ) + c 2 nj - c 2 2 n 2 . (47) 

We see here that as long as the momenta of the two com¬ 
ponents are not exactly parallel, i.e., ki is not exactly 
equal to k 2 , then |k| < |ki|, and there is always dynam¬ 
ical instability for the long wavelength excitations. 

Therefore, the spin current in Eq. (1461) is generally un¬ 
stable and not a super-flow. This instability originates 
from the interaction process described by V’oV’oV’iV’-i in 
the second quantized Hamiltonian. This energetically fa¬ 
vored process converts two particles in the to = 1,-1 
components, respectively, into two stationary particles in 
the to = 0 component. To suppress such a process and 
achieve a stable pure spin current, one can utilize the 
quadratic Zeeman effect. With the quadratic Zeeman ef¬ 
fect of negative coefficient, the Hamiltonian adopts an 
additional term A to 2 (A < 0 and m = 1,0,—1). This 
term does not change the energy of the m = 0 compo¬ 
nent, but lowers the energy of the other two components 
m = 1,-1. As a result, there arises a barrier for two 
atoms in the m = 1,-1 components scattering to the 
m = 0 component, and the scattering process is thus 
suppressed. 

The above intuitive argument can be made more rigor¬ 
ous and quantitative. Consider the case ki = k 2 . With 

the quadratic Zeeman term, the excitation spectrum for 
the to = 0 component changes to 


+ C 2 n — Aj —c\n 2 . (48) 

So as long as — A— Ii 2 |ki| 2 /2M > 0, long wavelength exci¬ 
tations will be stable for the m = 0 component. From the 
excitation of the m = 0 component, one can obtain a crit¬ 
ical relative velocity of the spin current, Uo = 2yJ—2X/M. 
There is another nonzero critical velocity v\ = 2 TIC 2 /M 
determined by the excitations of the to = 1,-1 compo¬ 
nents. The overall critical velocity of the system is the 
smaller one of vq and v\. Therefore, below the critical 
relative velocity v c = min{uo,Ui}, the pure spin current 
is stable and a super-flow. The experimental scheme to 
realize such a Zeeman effect will be discussed in subsec¬ 
tion [VC] 



B. Circular flow 

In the cylindrical geometry, we consider a pure spin 
current formed by two vortices with opposite circulation 
in the m = 1,-1 components. From similar arguments, 
one can expect that interaction will make such a cur¬ 
rent unstable. Inspired by the quadratic Zeeman effect 
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method above, we propose to use SOC to stabilize it. 
The SOC can be viewed as a momentum-dependent ef¬ 
fective magnetic field that exerts only on the m = 1,-1 
components. Therefore, it is possible that SOC lowers 
the energy of m = 1,-1 components, and consequently 
suppresses the interaction process leading to the insta¬ 
bility. 

The model of spin-1 BEC subject to Rashba SOC can 
be described by the following energy functional, 


£ [V’a] = dr \ 

•' l a= —1,0,1 


h 2 \Vlp a \ 2 C 0 2 C 2 2 


2 M 

+ 7 {SxPy - S y p x ) >, 


+ pV(r) + ^p 2 + -s- 


(49) 


(al) 


(bl) 


(cl) 


@ 


■2 0 2 
X 

( 82 ) 


0 


x 

(b2) 


® 


X 

(C2) 



where p is the density, V(r ) = \Mu> 2 (x 2 + y 2 ) is the 
trapping potential, and 7 is the strength of SOC. (• • •) 
is the expectation value taken with respect to the three 
component wave function %p = (ipi, tpo, %p_i) T . The SOC 
strength 7 defines a characteristic length a soc = h/Mj, 
and can be rescaled to be dimensionless with respect to 
the harmonic oscillator length ah = y/h/Mui. Then we 
characterize the strength of SOC with the dimensionless 
quantity k = ah/a soc = 7 The SOC of Rashba 
type here can be generated in various ways, which will 
be discussed in the next subsection. 

The above model can describe a spin-1 BEC of 23 Na 
confined in a two-dimensional harmonic trap. Assume 
the atom number is about 10 6 . Using the estimate of 
scattering lengths a 0 = 50as, 02 = 55as 74{, with clb 
being the Bohr radius, the ground state of spin -1 23 Na 
should be antiferromagnetic because Co > 0 ,C 2 > 0 [I. 7s]. 
Previous studies of spin -1 BEC with Rashba SOC mostly 
focus on the strong SOC regime, where the ground state 
is found to be the plane wave phase or the stripe phase, 
for ferromagnetic interaction and antiferromagnetic in¬ 
teraction, respectively [36j]. Here we are interested in 
the antiferromagnetic interaction case and the weak SOC 
regime (k -C 1 ), and calculate the ground state wave 
function of the energy functional with the method of 
imaginary time evolution. 

We find that when the SOC is weak (k <C 1), the 
ground state wave function has the form 


( Xi (r)e-*+ \ 

= Xo{r) , (50) 

V X-i (r)e» ) 


FIG. 8: (color online) Amplitudes (al,bl,cl) and phase an¬ 
gles (a2,b2,c2) of the three component wave function %p = 
(i/j 1 ,tpo,ip-i) T at the ground state of Hamiltonian (14911 for a 
BEC of 23 Na confined in a 2D harmonic trap. The particle 
number is 10®, the frequency of the trap is 2n x 42 Hz, and 
the dimensionless SOC strength is k = 0.04. The units of the 
X and Y axes are ah- 


of spin and angular momentum, the SOC term reads 


H 


SOC 


7 y/MHu 
2 




where S± is the ladder operator of spin, and a~\,(R) is the 
creation operator of the left (right) circular quanta [75$ . 
When the SOC is very weak (k <C 1), its effect can be 
accounted for in a perturbative way. From the ground 
state = | 0 , 0 ), the first order correction to the wave 
function for small 7 is given by 


vj/l 1 ) 


=2 ^(- S+a , +S - 4 )lo i0) 

= f (-| 1 ,- 1 ) + |- 1 , 1 », 


(52) 


where |m s ,ni 0 } denotes a state with spin quantum num¬ 
ber m s and orbital magnetic quantum number m 0 . One 
immediately sees that i/’i has angular momentum —h and 
ip -1 has angular momentum h. Besides, the amplitudes 
of both ipi and ip_ 1 are proportional to n. 

There exits a continuity equation for spin density and 
spin current, which is 


|(^S^)+V-J* =0. (53) 


with Xi(r) = —x_i(r) and all \i rea h The ground state 
is shown in Fig. [ 8 ] Such a ground state consists of an 
anti-vortex in the m = 1 component and a vortex in the 
m = —1 component. The m = 0 component does not 
carry angular momentum. Since \ip\\ = |-0_i|, the net 
mass current vanishes. 

The wave function in Eq. (1501) can be understood in 
the single particle level. In terms of the ladder operators 


The spin current density tensor Jf, (u, = x,y,z denotes 
the spin component) is defined as [76|, 177| 

=\ {^Spvip + c.c.} 

= \ | ^2 ‘•I’m ( S »)mn V ril^l + C.C. > , (54) 

I m,n,l J 
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where 

v n i = 1 (z X S ni ) , (55) 

and c.c. means the complex conjugate. The second part 
in is induced by the SOC. 

By the definition in Eq. (I54|l . the spin current density 
carried by the ground state (l50l) is 

J* = 7 sin 2 </)|V’i| 2 i + 7 (l^o| 2 +2|V’i| 2 sin 2 (f)y, 

= ~7 (I" 00 1 2 + 2|'0i| 2 cos 2 <f>) x - 7 sin 2 <()|^i \ 2 y, 

J- = p * 1 + V / 27lV ; i'0o|) </>• (56) 
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From both analytical and numerical results of the wave 
function, |-0i| <C |^o|i so J® roughly points in the y di¬ 
rection, while J® almost points in the —x direction. J® 
represents a flow whose amplitude has rotational symme¬ 
try. From the numerical results shown in Fig. [9] we see 
that J* is a counter-clockwise flow. The amplitudes of 
J* and J* are of the same order, both proportional to /c, 
while that of J®, proportional to k 2 , is much smaller. It 
is evident that the state in Eq. m carries no mass cur¬ 
rent and only pure spin current. Since the spin current 
is in the ground state, it must be stable. In this way, we 
have realized a superfluid of pure spin current, or a pure 
spin super-current. 

C. Experimental schemes 


2 

1 

> 0 
-1 
-2 


-2 

-1 

0 

1 

2 



X 





(c) 




2 


In this subsection, we propose the experimental 
schemes to generate and detect the pure spin currents 
discussed in subsections IV Al and IV Bl 

The planar pure spin current can be easily generated. 
By applying a magnetic field gradient, the two compo¬ 
nents m = 1,-1 will be accelerated in opposite direc¬ 
tions and _a pure spin current is generated as done in 
Refs. [H,j63]. To stabilize this spin current, one needs to 
generate the quadratic Zeeman effect. We apply an oscil¬ 
lating magnetic field B sin uit with the frequency co being 
much larger than the characteristic frequency of the con¬ 
densate, e.g., the chemical potential g. The time averag¬ 
ing removes the linear Zeeman effect; only the quadratic 
Zeeman effect remains. The coefficient of the quadratic 
Zeeman effect from the second-order perturbation theory 
is given by A = ( g^B ) 2 /AEhf, where g is the Lande g- 
factor of the atom, /ie is the Bohr magneton, and AEfif is 
the hyperfine energy splitting [7^ ] . For the F = 2 mani¬ 
fold of 87 Rb, A E hi < 0, so the coefficient of the quadratic 
Zeeman effect is negative. 

The circular flow in subsection IV Bl may find prospec¬ 
tive realizations in two different systems: cold atoms and 
exciton BEC. In cold atoms, we consider a system con¬ 
sisting of a BEC of 23 Na confined in a pancake trap, 
where the confinement in the 2 direction is so tight that 
one can treat the system effectively as two dimensional. 
The SOC can be induced by two different methods. One 
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-1 
-2 


- 2-1 01 2 

X 

FIG. 9: (color online) Distribution of the spin current densi¬ 
ties JJ (a), Jy (b) and J| (c) of the ground state shown in 
Fig. [8] The length of the arrows represents the strength of 
the spin current. The arrow length of different subfigures is 
not to scale, k = 0.04. The units of the X and Y axes are ah- 
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is by the exertion of a strong external electric field E in 
the z direction. Due to the relativistic effect, the mag¬ 
netic moment of the atom will experience a weak SOC, 
where the strength 7 = <7^3 |E| /Me 2 . Here M is the 
atomic mass and c is the speed of light. For weak SOC 
(small 7 ), the fraction of atoms in the m = 1,-1 com¬ 
ponents is proportional to 7 . For an experimentally ob¬ 
servable fraction of atoms, e.g., 0 . 1 % of 10 6 atoms, using 














































14 


the typical parameters of 23 Na BEC, the estimated elec¬ 
tric field is of the same order of magnitude as the vacuum 
breakdown field. For atoms with smaller mass or larger 
magnetic moment, the required electric field can be low¬ 
ered. Another method of realizing SOC is to exploit the 
atom laser interaction, where strong SOC can be created 
in principle [7f|. In exciton BEC systems, as the effective 
mass of exciton is much smaller than that of atom, the 
required electric field is four to five orders of magnitude 
smaller, which is quite feasible in experiments [80l483j . 

The vortex and anti-vortex in the m = 1,-1 compo¬ 
nents can be detected by the method of time of flight. 
First one can split the three spin components with the 
Stern-Gerlach effect. The appearance of vortex or anti¬ 
vortex in the m = 1,-1 components is signaled by a 
ring structure in the time of flight image. After a suffi¬ 
ciently long time of expansion, the ring structure should 
be clearly visible. 

VI. SUMMARY 

In summary, we have studied the superfluidity of three 
kinds of unconventional superfluids, which show distinct 
features from a uniform spinless superfluid. The peri¬ 


odic superfluid may suffer a new type of instability, the 
dynamical instability, absent in homogeneous case; the 
spin-orbit coupled superfluid has a critical velocity de¬ 
pendent on the reference frame, a new phenomenon com¬ 
pared with all previous Galilean invariant superfluids; the 
superfluid of a pure spin current, though scarcely stable 
in previous studies, can be stabilized by the quadratic 
Zeeman effect and SOC in planar and circular geometry, 
respectively. These new superfluids significantly enrich 
the physics of bosonic superfluids. 

With the rapid advances in cold atom physics and 
other fields, the family of superfluids is expanding with 
the addition of more and more novel superfluids. Pre¬ 
vious study has greatly deepened and enriched our un¬ 
derstanding of superfluidity, but we believe physics more 
exciting and beyond the scope of our current understand¬ 
ing remains to be discovered in the future. 
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